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ABSTRACT 

We present the first detection of neutral atomic oxygen ( 3 Pi- 3 P 2 fine structure 
line at ~63 /zm) toward the Horsehead photodissociation region (PDR). The cloud 
has been mapped with the Spitzer Space Telescope at far-IR (FIR) wavelengths 
using the Multiband Imaging Photometer for Spitzer (MIPS) in the spectral 
energy distribution (SED) mode. The [Ol]63/im line peaks at the illuminated 
edge of the cloud at y4y~0.1-0.5 (inward the gas becomes too cold and outward 
the gas density drops). The luminosity carried by the [Ol]63/im line represents 
a significant fraction of the total FIR dust luminosity {Iqs/Ifir— 4x10 -3 ). We 



lr rhis work is based on observations made with the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
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analyze the dust continuum emission and the nonlocal O I excitation and ra- 
diative transfer in detail. The observations are reproduced with a gas density 
of rin- 10 4 cm -3 and gas and dust temperatures of Tj.~100 K and K. We 

conclude that the determination of the O I 3 Pj level populations and emergent 
line intensities at such "low" densities is a complex non-LTE problem. FIR ra- 
diative pumping, [Ol]63/im subthermal emission, [Ol]145//m suprathermal and 
even maser emission can occur and decrease the resulting [Ol] 63/145 intensity 
ratio. The Herschel Space Observatory, observing from ~55 to 672 /zm, will allow 
us to exploit the diagnostic power of FIR fine structure lines with unprecedented 
resolution and sensitivity. 

Subject headings: astrochemistry — infrared: ISM — ISM: atoms — ISM: clouds 
— line: formation — radiative transfer 



1. Introduction 

Oxygen is the most abundant heavy element in the Universe. Due to its relatively 
high ionization potential, far-ultraviolet (FUV) photons (6 eV<hu<13.6 eV) do not ionize 
it and thus a large fraction of atomic oxygen can be expected in the photodissociated gas 
interface between the ionized and the molecular gas before chemical species such as CO, 
H2O or OH reach high abundances. The absorption of FUV photons by dust grains leads 
to the ejection of photoelectrons that heat the gas (e.g., Draine 1978) as well as to intense 
FIR dust continuum emission. The efficiency of this essential process can be estimated by 
comparing the luminosity of the main gas cooling lines (e.g., the FIR fine structure lines of 
Cll and Ol) with the luminosity of the dust emission. Indeed, FIR lines from abundant 
gas reservoirs are key players in the energy balance of the regions they trace (from FUV 
illuminated interstellar clouds and protoplanetary disk surfaces to extragalactic starbursts). 

The Ol( 3 Pj) fine structure triplet system is formed by the 3 P 2 (ground-state), the 
3 Pi (A£i 2 /£;=228K) and the 3 P (AE 01 /k=99K) levels. Radiative transitions between 
the 3 Pi- 3 P 2 and 3 Po- 3 Pi states lead to the widespread [Ol]63 and 145 /im line emission 
respectively. Observations and models of such lines are crucial diagnostics of a broad range 
of galactic and extragalactic environments (see Hollenbach & Tielens 1999 for a review). 

In this Letter, we present the first detection of the [O l]63 /im line toward the Horsehead 
PDR (<i~400pc, edge-on), one of the most famous FUV-induced chemical laboratories in the 
sky. We analyze the FIR line and continuum emission without strong simplifying assumptions 
in the excitation and radiative transfer treatment. 



- 3- 



2. Observations and Data reduction 

The Horsehead nebula was observed with the Multiband Imaging Photometer for Spitzer 
(MIPS) as part of our "SPECPDR" program (ID:20281; Joblin et al. 2005) on 2006 April 3 
using the spectral energy distribution mode (hereafter MIPS/SED). MIPS/SED provided 
the capability of obtaining low-resolution spectra (R = 15 — 25) in the FIR range (A = 
55 — 97/im) using a slit spectrometer and the MIPS 70 Ge:Ga detector array. The slit 
size has two detector pixels in width (~20") and 12 pixels in length (~2'). A reflective 
grating provides a spectral dispersion of 1.7/im pixel -1 . We used the "spectral mapping" 
mode and covered the Horsehead with 14 pointings shifted by half the slit width in the 
direction perpendicular to the slit length. The exposure time was 10 s per pointing. The 
"chop distance" for the OFF positions was set to +3' for each pointing. Both ON- and 
OFF-source areas are shown in Fig.[]}i. 

We subtract the average OFF-source spectrum to every ON-source spectra. We first 
retrieved the BCD data (ver. sl6.1.0) that we proceed with the mosaic_sed.pl script of 
MOPEX0 to obtain the post-BCD level data (in intensity units of MJysr We apply 
a flat correction to the spatial dimension on the detector (i.e. detector columns relative 
gain correction) at the BCD level. We compute the flat response of the detector columns 
using the OFF positions. We have developed a pipeline that builds spectral cubes (two 
spatial dimensions and one spectral dimension) from the Post-BCD data. We work in the 
assumption that no spatial information is available in the slit width. Thus, the N\ rows of 
the detector can be seen as images of the sky through the slit for A^a sampled wavelengths 
between 55 and 97 fim. For each wavelength, this image of the slit is reprojected into a 
N a x Ng spatial grid of the sky. The final cube then has N a x N$ x N\ pixels. The pixel 
size of the spatial grid is half of each detector pixel field-of-view, PFOV (9.8" /2 = 4.9"). 

Nevertheless, due to dispersion, the light that falls on a pixel comes from the entire slit 
width, that is two times the pixel size, so that the PFOV is rather 9.8" x 19.6" . This leads to 
an overestimation of the intensity of a factor 2 that we take into account here. For extended 
sources, one also has to correct for the effect of the aperture used for the photometric cali- 
bration with point sources. We then divide the intensity by the 5 pixels aperture correction 
function (see MIPS Data Handbook). The photometry of the resulting spectral cubes (for 
all sources in our program) were checked by cross- calibrating them with MIPS 70 data when 
available. It appears that the MIPS/SED gain is 30% higher than the MIPS 70 one. On the 
basis that the photometer calibration is more reliable, we apply this 30% correction to our 
MIPS/SED data and use the MIPS 70 absolute calibration uncertainty of 7%. 



http://ssc.spitzer.caltech.edu/postbcd/mopex.html 
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3. Observational Results 

FigurefTJx shows the ON- and OFF-source areas of our Spitzer observations plotted 
over the polycyclic aromatic hydrocarbon (PAH) mid-IR emission toward the Horsehead 
(Abergel et al. 2003). The PAH band emission traces the FUV illuminated edge of the 
cloud (Compiegne et al. 2007, 2008; Goicoechea et al. 2009). MIPS spectra averaged over 
the ON and OFF areas are shown as green and blue histograms (Fig.QJ). The ON-OFF 
subtracted spectrum is shown as a black histogram. The [Ol]63/im line is not present in 
the OFF spectrum, and therefore, the extended line emission can fully be attributed to the 
cloud. 

In order to have an approximate view of the spatial location of the oxygen emission, 
Fig.[Tb displays the [Ol]63/im line intensity measured in each pixel of the data cube (gray 
scale). The ortho-H.2 v—0-0 S(l) pure rotational line emission at ~17/im, delineating the 
edge of the PDR, is shown as white contours. Overall, the [O i] emission peaks in the warm 
cloud layers (where H2 and PAH emission are also bright) directly exposed to the FUV 
radiation from cxOri 09. 5V star. Although the excitation requirements of the ortho-H.2 v=0- 
S(l) (Ejk ~1015K) and [Oi]63 /mi (E u /k ~228K) lines are different, the apparent shift 
between both line emission peaks in Fig. lc may be marginal since MIPS/SED astrometry is 
not better than ~5". Besides, cloud inclination variations (relative to the line of sight) from 
north to south, and thus differential line opacity effects, likely explain the brighter [O l]63 /im 
emission observed toward the north of the filament. 

The [O l]63 /mi line intensity toward the emission peak, as measured over the ~15" full 
width at half maximum (FWHM) angular resolution of MIPS at 63 /im (circle on Fig. [lb), is 
J 63 =(1.04±0.14)xl0~ 7 Wm^sr 1 (~4 times brighter than the H 2 v—1-0 S(l) fluorescent 
emission at 2.12 /mi; see Habart et al. 2005). The resulting MIPS spectrum is shown in Fig.© 
as "PEAK". The [Ol]63/mi line is detected at a 10<r level. The FIR thermal emission of 
dust raises from ~155 MJy sr _1 (at ~55 /im) to ~325 MJy sr _1 (at ~85 /im). The continuum 
peaks at longer wavelengths and thus the temperature of the grains emitting in the FIR (T<j) 
cannot be high (T d <35K), lower than the expected gas temperature (T k ). Apart from this 
[Ol]63/mi bright emission toward the cloud edge, the MIPS data show a fainter level of 
widespread emission (~4 times weaker) likely arising from the more diffuse cloud envelope, 
or halo, illuminated by the ambient FUV radiation field. 
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4. Discussion 

4.1. Physical conditions in the Horsehead PDR 

In order to guide our interpretation of the [Ol]63/im line toward the PEAK, Fig. [2] 
(left) shows the output of a photochemical model (Le Petit et al. 2006; Goicoechea & Le 
Bourlot 2007) adapted to the steep density gradient, nn~r -3 (Habart et al. 2005, Goicoechea 
et al. 2009) and FUV illumination conditions in the Horsehead edge (Go— 100, where Go is 
the FUV flux in Habing units of the local interstellar value). We adopt an elemental oxygen 
abundance of 3xl0 -4 relative to H nuclei (Savage & Sembach 1996). Taking into account 
that the gas temperature rapidly decreases from the low density PDR (T^^IOOK; see also 
Pety et al. 2005; Habart et al. 2009) to the denser and shielded inner layers (T fc <20K; Pety 
et al. 2007), we predict that the [Ol]63/im emission peak should occur at the surface of the 
molecular cloud between Ay— 0.1 and 0.5. Due to the relatively low FUV radiation field 
impinging the cloud, dust grains are not expected to attain large temperatures (T^<30K). 
However, H 2 photodissociation is still important in the cloud edge and a significant fraction 
of hydrogen is predicted to reside in atomic form in the same layers where [Ol]63/im is 
bright. At extinction depths greater than Ay~1.5, carbon monoxide starts to be protected 
against photodissociation and the gas temperature falls sharply. The formation of CO (and 
other O-rich gas and ice species) decreases the abundance of atomic oxygen inward. 

4.2. Nonlocal and Non-LTE Modeling 

Observations of the [O l]63 /im line and FIR continuum emission constrain the mean gas 
density and both the dust and gas temperatures toward the PDR. In this work, we use a 
nonlocal and non-LTE code that treats both the line and continuum radiative transfer and 
includes the radiative pumping of O I levels by absorption of FIR dust photons (see the 
Appendix in Goicoechea et al. 2006). Almost all studies solving the excitation of atomic 
fine structure levels simplify the problem by assuming a local treatment of the excitation 
calculation: LVG or single-zone escape probability methods (e.g., Rollig et al. 2007), where 
the source function (the line emissivity/absorption ratio) is assumed to be constant, and the 
radiative coupling between different cloud positions is neglected. The effects of dust opacity 
are also generally omitted (Gonzalez Garcia et al. 2008). Their applications to the inter- 
pretation of [O l]63,145/zm lines have been widely used by a number of authors in the past 
(e.g., Tielens & Hollenbach 1985; van Dishoeck & Black 1986; Spaans et al. 1994; Poglitsch 
et al. 1995; Kaufman et al. 1999; Vastel et al. 2001; Liseau et al. 1999, 2006). Nevertheless, 
Elitzur & Asensio-Ramos (2006) recently compared the exact solutions provided by their 
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nonlocal method with those given by the local approximation and concluded that, for gas 
densities around 10 3 -10 4 cm -3 , the error in the predicted [Oi] 63/145 line intensity ratio can 
be significant (i. e., comparable with the variation range of the ratio for different temperatures 
and densities) . Indeed, these are the densities expected in the outer layers of the Horsehead 
nebula, i.e., below the critical densities of Oi fine structure transitions (n cr ^5x 10 5 cm~ 3 
for O I ( 3 Pi- 3 P 2 )-Il2 collisiona@ de-excitation). Under these conditions, the determination 
of the 3 Pj level populations is a non-trivial problem and radiative and opacity effects (FIR 
pumping and line-trapping) can play a role. Besides, the [Ol]63/im line becomes optically 
thick (t63 > 1) for low column densities, N(0 1)>10 17 cm -2 , and the [Oi]145/xm line can be 
suprathermal (T ex >T k ) or even inverted (T ex .<0) and the emission amplified (r 145 <0) by the 
maser effect (see also Liseau et al. 2006; Elitzur & Asensio-Ramos 2006). 

4.3. FIR dust continuum emission 

In order to constrain the dust temperature, the FIR dust luminosity and obtain an ap- 
proximated value for the gas density in the [O l]63 /zm emitting region, we have first modeled 
the observed SED toward the PEAK position. We explicitly solve the continuum radia- 
tive transfer problem by taking a dust opacity r^(A) proportional to k 10 o(100/A)^, where 
Kioo is the dust opacity at 100 /xm (per gas+dust mass column density) and (3 is the dust 
spectral index. Assuming a line-of-sight depth of 0.1 pc (uncertain within a factor of ~2, 
see constraints by Habart et al. 2005), the best models of the continuum shape and ab- 
solute intensity level are obtained for Td—30 K (in agreement with the PDR model predic- 
tion) and nH=Ti(H)+2n(H2)— 10 4 cm" 3 (with /3~1.4 and /tioo— 0-25 cm 2 g -1 ; e.g., Li & Draine 
2001). Integrating the resulting continuum intensity distribution from 40 to 500 fim we ob- 
tain Jfir=2.5x10 -5 Wm~ 2 sr" 1 . The dust model is shown in Fig.© (red dashed curve). We 
determine that the line-of-sight extinction depth toward the PEAK is Ay— 1-5, consistent 
with the expected values of a cloud edge. The FIR dust continuum toward the region of 
bright [Oi], H 2 , and PAH emission is optically thin, with Td(63)~2xl0 -3 . 



2 The collisional rate coefficients (7^ ) used in this work are: 1 collisions with ortho /para-H2 from Jaquet 
et al. (1992): note that for low energies (Tfc<300-400K), the 702 and 712 rates are slightly higher (e.g., by 
<10% at 100 K) when para-H2 is the perturber. Otherwise, rate coefficients are slightly higher for ortho-H2 
collisions; Oi-H collisions from Abrahamsson et al. (2007): at 100 K these excitation rate coefficients are a 
factor of ^4 larger than those by Launay & Roueff (1977); O i-He collisions from Monteiro & Flower (1987); 
Oi-e" collisions from Bell et al. (1998) and Oi-p + collisions from Chambaud et al. (1980). 
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4.4. O I excitation and radiative transfer 

After modeling the FIR radiation field in the PDR we now solve the O I 3 Pj populations 
adding radiative pumping by FIR dust photons. We include collisional rate coefficients 2 
(cm 3 s -1 ) of atomic oxygen with ortho- and para-H.2, H, He, electrons, and protons. H 2 and 
H collisions dominate, but it is important to include their rates (s _1 ) with the most realistic 
H and H2 density contributions. Note that, for example, 3 Pi — 3 Pq collisional excitation can 
be neglected (to first order) if only O1-H2 collisions apply (Monteiro & Flower 1987), but it 
does contribute if collisions with H are included (Abrahamsson et al. 2007). We adopt an 
ortho-to- para-H.2 ratio 2 (OTP) of 3, although, at the and nn considered in this work, our 
results do not change by taking an OTP of 1 ([Ol] lines become ~2% brighter). 

Figure [2] (right) displays a grid of nonlocal models showing the emergent Iqs inten- 
sity as a function of N(0 1) column density. The lower panel shows the predicted Iqs for 
7^=100 K (i.e., roughly the gas temperature in the outermost cloud layers) and varying 
H 2 densities. Guided by our PDR model (Figure [2] (left)), and in order to modify only 
one physical parameter in the grid, we take a fixed n(H) density of 10 3 cm -3 and a high 
ionization fraction: n(e~)/nn=10~ A (see Goicoechea et al. 2009). Our calculations in- 
clude thermal, turbulent, and opacity line broadening with a turbulent velocity dispersion 
of a=0.225 km s _1 (FWHM= 2.355 x a). This value is inferred from heterodyne interfero- 
metric observations of the HCO radical toward the PDR (Gerin et al. 2009). The shaded 
regions show the observed [O l]63 /mi line intensity toward the PEAK position. For a line-of- 
sight depth of 0.1 pc and taking Ol/H~3xl0~ 4 in the PDR (most oxygen in atomic form), 
our observations are reproduced with n(H 2 )~5xl0 3 cm -3 (consistent with the dust model) 
and N(0 l)~8.5xl0 17 cm -2 (best model shown in Figure 1 (b)). The [Ol]63/mi opacity at 
line center is r 63 ~6 (due to saturation the width of the line profile is opacity broadened, 
FWHM~1.3kms _1 , i.e., broader than the linewidth of an optically thin line affected only 
by turbulent and thermal dispersion). The predicted [Ol]63/145 line intensity ratio is ~15. 
The required H 2 density has to be multiplied (divided) by ~4 if the gas temperature is a 
factor of 2 lower (higher); line opacities and [Ol] 63/145 intensity ratios change accordingly. 
Future observations of both [0 1] 63, 145 /im lines will help to better constrain the density and 
the atomic oxygen abundance gradient. 

For Tfc=100K, the [Ol]63/nn line emission is subthermal (T ex <Tk) which means that 
the 3 Pi level (n x ) is much less populated than in LTE (i.e., ^ <C |e~^§ In fact we 

find that the 3 P - 3 Pi transition at 145 /im is suprathermal (with ^ ~ | > |e~^> ~|). This 
effect is produced by the larger radiative life time of the upper 3 P level compared to the 
lower 3 Pi level (by a factor of ~5). Even more, the lack of significant 3 Po- 3 Pi collisional 
de-excitation for O l-H 2 collisions (neglecting those with H) can invert the level populations 
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(m > I) leading to a [Ol]145/im maser line emission. High temperatures (T fc >100K for 
a pure H 2 gas and Tfc>300 K for a pure H gas) together with relatively low r 63 opacities 
(or nH<10 5 cm~ 3 ) are required to obtain the level inversion, otherwise line-trapping in the 
[Ol]63/im line increases the upper 3 Pi level population again. At such "low" gas densities 
(bh < n CT ) the local treatment of O I excitation overestimates (underestimates) the ri\ (n ) 
level populations, leading to stronger (weaker) [O l]63 /itn (145 /zm) emission as long as t 63 >1. 
The error made increases with iV(Ol) (see Elitzur & Asensio-Ramos 2006). 

The FIR continuum toward the Horsehead is optically thin and radiative pumping does 
not play a significant role in the Ol excitation (it reduces I$3 by ~1 %). However, in more 
extreme sources (with opaque continuum emission), FIR pumping can reduce the T ex of the 
Ol 3 Pi- 3 P 2 transition (T ex ^ T d ) and increase the n /ni level population ratio. For the 
physical conditions treated here, the [O i] 63/ 145 ratio decreases by a factor of ~2 if the dust 
emission becomes optically thick at 63 /zm. Overall, neglecting FIR pumping and using local 
approximations for the Ol excitation overestimates the [Ol]63/145 ratio. 

4.5. Cooling rate and photoelectric heating efficiency 

Once the non-LTE level populations have been accurately computed, one can determine 
the appropriate [O l]63 /itn cooling rate per volume unit, A 63 =/iz/i 2 (n 2 C 2 i-ni Ci 2 ) (where 
are the collisional excitation and de-excitation rates in s _1 ). This relation takes into account 
that only line photons following collisional excitation and escaping the cloud are responsible 
of the gas cooling. We compute A 63 ~5xl0~ 21 ergs -1 cm -3 toward the PEAK position. 

Our observations and models are consistent with the [Ol]63/im line arising from the 
outer layers (A v ^0.1 and 0.5) of the Horsehead PDR (G — 100 and %- 10 4 cm~ 3 ) where the 
gas is warm (~100K) and not all available hydrogen has converted into H 2 yet. Together 
with available [C n] 158 /iin observations (Ii 58 — (1— 3) x 10~ 7 Wm -2 sr -1 , from KAO at ~55" 
resolution; Zhou et al. 1993), we find that the luminosity of the [Cll]158/im + [Ol]63/im 
lines represents ~l%-2% of the total FIR dust luminosity. This fraction corresponds to a 
high efficiency of the photoelectric heating mechanism. In fact, due to the relatively low FUV 
field but moderate gas density in the PDR, the Horsehead is a good source to investigate 
the grain populations that efficiently transform the incident FUV field into gas heating. 

Herschel will allow us to map almost all relevant FIR cooling lines at much higher 
spectral and angular resolutions. Together with appropriate models of their excitation, it 
will then be possible to study the energy balance of interstellar clouds in greater detail. 
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Fig. 1. — (a) ISOCAM (5-8.5 /im) image of the Horsehead (dominated by PAH emission; 
Abergel et al. 2003). The ON- and OFF-source areas of our Spitzer /MIPS observations 
are shown, (b) Different spectra: mean ON-source (green histogram), OFF-source (blue 
histogram), and ON-OFF (black histogram). The MIPS spectrum toward the PDR (PEAK 
position) is shown on top (black continuous curve) together with our dust and [O i] radia- 
tive transfer model convolved with MIPS/SED spectral resolution, (c) [Oi]63/xm inten- 
sity map. The scale is shown in the look-up table (in Wm^sr 1 ). The white contours 
represent the ortho-H.2 v=0-0 S(l) line intensity at ~i7/im (0.9 x, 1.5 x, 2.5 x, 3.5 x, and 
4.0xl0 -8 Win -2 sr -1 ; Habart et al. 2009). The circle represents the ~15"-FWHM resolution 
toward the PEAK. 
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Fig. 2. — Left: PDR model adapted to the Horsehead (Go=100 and density gradient from 
Goicoechea et al. 2009). The upper panel shows the resulting H 2 and H density profiles, 
together with the Tk and Td temperature gradients. The lower panel shows the CO/C/C + 
transition and the predicted abundance profile of atomic oxygen. Right: output of a grid 
of nonlocal, non-LTE radiative transfer models. Panels show the expected [Ol]63/im line 
intensity (in Wm _2 sr _1 ) as a function of iV(Oi) for two gas temperatures (Tfc=50 and 
100 K) and different gas densities: n(H)=10 3 cm -3 fixed and varying n(H 2 ). The shaded 
regions show the [Ol]63yum line intensity observed by Spitzer /MIPS toward the PEAK. 



